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Abstract 

We demonstrate the classical stability of the weak/Planck hierarchy within 
the Randall-Sundrum scenario, incorporating the Goldberger-Wise mecha- 
nism and higher-derivative interactions in a systematic perturbative expan- 
sion. Such higher-derivative interactions are expected if the RS model is 
the low-energy description of some more fundamental theory. Generically, 
higher derivatives lead to ill-defined singularities in the vicinity of effective 
field theory branes. These are carefully treated by the methods of classical 
renormalization. 
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I. INTRODUCTION 



Theories with extra dimensions provide new ways of explaining the weak/Planck hierar- 
chy. The original proposal for doing so appeared in Ref. An alternative proposal is the 
Randall-Sundrum (RSI) scenario [0, where the hierarchy is set by a relative warp factor, 
^-kirrc^ between a "visible" brane, to which the Standard Model is confined, and a "hidden" 
brane where 4D gravity is highly localized by the RS2 mechanism 0. Here, /c is a fundamen- 
tal scale determined by the 5D cosmological constant and Vc is the compactification radius. 
The Goldberger-Wise mechanism Q provides a simple and natural means of stabilizing the 
radius at kvc ~ 0{l/e), by introducing a bulk scalar field with 5D mass-squared of order 
e in fundamental units. The large observed weak/Planck hierarchy, e~'^'^^'= ~ 10"^^, is then 
generated from a modest fundamental hierarchy, e ~ 1/10. 

Since the RSI field theory, including general relativity, is quantum- mechanically non- 
renormalizable, the model must be considered to be an effective description of a more fun- 
damental theory. Refs. have discussed string theory embeddings of the RSI mechanism. 
In any such embedding, higher-derivative interactions {a' corrections in string theory) are 
expected to appear in the effective field theory after integrating out very massive physics. It 
is, therefore, important to demonstrate that RSI and the Goldberger-Wise mechanism are 
stable under the addition of such higher-derivative terms. In this paper, we will show that 
this is indeed the case within the systematic framework of classical effective field theory. 

While short- distance quantum effects can be parameterized and studied within a local 
derivative expansion, it is also important to demonstrate stability in the presence of gen- 
uine, long-distance quantum effects. These can also be studied within effective field theory. 
Refs. [^] have examined such effects at one loop. We hope to give a more complete treatment 
of quantum effects in future work. 

Recently, a dual picture of the RS scenario has been developed [0], based on the AdS/CFT 
correspondence p[, in which the extra-dimensional dynamics is replaced by a strongly cou- 
pled conformal field theory. While this duality is compelling and powerful, we will not make 
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use of it in this paper as some aspects remain unproven. 

Our strategy is to first set up a systematic perturbative expansion for the classical ef- 
fective field theory in which to study higher derivatives. While the Goldberger-Wise e is an 
obvious small parameter of the scenario, we cannot perturbatively expand in it since the 
hierarchy is set by e~^/^, which vanishes to all orders in e. Instead, we choose both the 
bulk curvature and the brane tadpole couplings which give the Goldberger-Wise scalar a 
non-trivial profile in the extra dimension to provide our formal expansion parameter, A. In 
Section II, we re-derive the Goldberger-Wise mechanism in the absence of higher derivatives. 
We note that there is an elegant, exactly soluble version of the Goldberger-Wise mechanism 
P], but our perturbative treatment will be more convenient when higher-derivative terms 



are added. A discussion of the Goldberger-Wise mechanism related to ours is Ref. [10 



In Section III, we discuss how higher-derivative terms are constrained by symmetries. In 
Section IV we discuss the apparent incompatibility of the derivative expansion, normally 
valid at long distances, with the presence of "thin" or 5-function branes. We show how the 
ill-defined singularities that arise in the equations of motion can be eliminated by classical 
renormalization. In Section V we demonstrate the stability of the RS and Goldberger-Wise 
mechanisms when higher-derivative perturbations are included. The central technical con- 
cern is that terms in our perturbative expansion take the form A"/„(e) and it is important for 
a controlled expansion that the small parameter A is not overwhelmed by possible large terms 
in fn, such as e^^" or l/e"^. This is carefully checked. Section VI provides our conclusions. 



II. THE GOLDBERGER-WISE MECHANISM IN PERTURBATION THEORY 

A. The model 

RSI has a single extra dimension which is an interval, realized as an orbifold, S^/Z2. We 
will begin by using a conventional angular coordinate, — vr < < vr, for the S^, where the 
orbifold symmetry acts by ^ —0. We will always describe fields within the fundamental 
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domain < < tt. Their extension to general (p is then determined by the orbifold symmetry 
and periodicity on S^. 

The Goldberger-Wise mechanism will first be implemented within a theory given by 

S = Shuik + Syis + Shid, (2.1) 

where 

Sbuik = J d'x £^ dcj>VG (-^R + ^{Dxf + 3A;2 - ^ek^^^ (2.2) 

S^is = -M^k J d^x^ (^p^ + Kx + \l^vX'^ (2-3) 

Shid = -M^/c y d^x^ {^ph + AftX + ^At/^X") ' (2-4) 

with 

9lv{x) = G^^{x,(t)^ tt) 

9lA^) = G^,{x,ct> = Q). (2.5) 

Note that we have chosen a normalization such that x is dimensionless. We assume that 
there are no extremely large hierarchies among the couplings of the model. However, we do 
take 

k<M 

\K,h\ < e < \pv,h\ ~ Pv,h ~ 1. (2.6) 

We will restrict our attention to classical solutions which admit a four- dimensional 
Poincare invariance. Such configurations satisfy the ansatz 

ds^ = e^^^'t'^ri^^dx^dx" - rld(f>^ 
X = X(0), (2.7) 

where Tc is the (constant in this ansatz) "radius" . In solving the equations of motion it is 
convenient to work with a re-scaled dimensionless extra-dimensional coordinate, 



y = kr,(t), (2.8) 

so that the infinitesimal distance in the extra-dimension is dy/k. The Poincare ansatz then 
reads, 

ds^ = e^^^y^r]f,^dx^'dx'' - dyyP 
x = xiy)- (2.9) 

The equations of motion subject to this ansatz are, 

x" + 4A'x' -ex = {K + fivx)S{y - yc) + + iJ^hx)Ky) (2- lo) 

2 2 1 

^' + -^{xf = -3 (P^ + ^vX + ^f^vX^)S{y - yc) 

-l{ph + Kx + ^f^hX^)S{y) (2.11) 
{Ar = l-lex' + l{xr, (2.12) 

where 

?/c = kVcTT. (2.13) 

The brane-locahzed scalar tadpoles provide us with our formal small expansion param- 
eters, Xy^h ~ We expand our solution as a perturbation series in A, 

n 

X = EXn, (2.14) 

n 

where the subscript n denotes the term of order A" in the series. Notice that at order 
A° we have x = and the (truncated) AdS^ gravity solution obtained in RSI, Aq = —y, 
with Tc arbitrary. Thus, in fluctuations away from 4D Poincare invariance, becomes a 
"radion" modulus at zeroth order. We will demonstrate the Goldberger-Wise mechanism 
for stabilizing the radius in higher orders of perturbation theory. 

The strategy for solving the equations of motion is as follows. We solve for A' in terms 
of X using ( P.12|) , 
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(2.15) 



and eliminate A' from (|2.10|) to obtain an equation purely for x, 



X" - 1 - \ex^ + \{xT -ex = {K + ^^vX)Ky - + {K + (2.16) 



We will solve ( p.l6|) to any desired order in A. We then integrate (|2.15|) to solve for Aiy), 



subject to the canonical gauge choice ^(0) = 0. Equation ( p.ll|) will then be automatically 



solved away from the branes as a consequence of 5D general covariance.Q Finally, we satisfy 



the two brane junction conditions of ( 2.11 ) by fine-tuning the hidden brane tension param- 
eter ph (equivalent to fine-tuning the 4D cosmological constant to zero), and adjusting the 
compactification radius, Tc (or equivalently, Uc) to its stable vacuum value. 

B. Perturbation theory 
Here, we will show self-consistently that the solution of the equations of motion satisfies, 

A{y) ^ -1, (2.17) 

so that the weak/Planck hierarchy determined by the RSI mechanism is set by e~'~^^^f'^\ 
Therefore, while we need e to be somewhat small, we cannot work perturbatively in e. (As 
stated above we will be working strictly perturbatively in A.) However, we will simplify 
our analysis by dropping subleading terms in e~'^'^^/'^\ Note that while A is formally small, 
combinations such as A(9(l/e'") or A(9(e~''^/^) may be large when e ^ 1. It is, therefore, 
crucial that such combinations do not appear in the perturbative series in order for the 
perturbative expansion to be under control. We show that this danger does not eventuate 
by a careful analysis. 



^One can easily check that up to 5-function terms, ( 2.1C| ) and (|2.12|) imply ( 2.1l|) . 
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1. Stabilization at second order 



Since X vanishes at zeroth order, we can perturbatively expand the square- root in ( p. 16 ). 
At first order, 

Xi - 4x'i - exi = (A, + iJ,xi)S{y - Vc) + {Xh + i^hXi)S{y)- (2.I8) 

With orbifold boundary conditions the solution is 

^^ = c,e^+iy-yc) ^C2e^-y^ (2.19) 



with 



and where 



A± = 2 ± yiTe. (2.22) 



Note that for small e, 



A+ ^A, A_ ^ -e/4. (2.23) 

In the expressions for the coefficients we have neglected subleading powers of e~^+^'', since 
we will demonstrate (p.l7|) . 



At first order in A the scalar field does not stabilize the radion because A receives no 
correction at this order. There are still two fine tunings needed to satisfy the equations of 
motion as in the original version of RSI without stabilization. It is only at second order, 
where the first order scalar profile back-reacts on the metric, that the compactification radius 
is fixed. This back-reaction was also discussed in Ref. |]TU|. At second order X2 = and A'2 



is given by ( |^.12| ) expanded to second order in %. The junction conditions of ( |2.11| ) at this 
order read, 
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(-1 + ^%(0)' - ^XM') = -I (p. + A,xi(0) + ^P.Xi(O)^) (2.24) 

(-1 + ^^XiiVcf - ^x'i(yc)') =l{p- + ^vXiiVc) + ^fivXiiVc?^ ■ (2.25) 

Note that parametrically in A, the only way the visible junction condition can be solved 
is if 

= 3 + p, ~ 0(A2). (2.26) 

Such a condition will not reappear at higher orders in A. Once we grant that pv is somewhere 
in this (9(A^)-sized window about —3, the visible junction condition can be satisfied, not by 
fine tuning of couplings, but by solving for the stable vacuum value of the dynamical radius, 

V.^'-^. (2.27) 

where 

A+(2A_-/Xfe)(2A_-2-/i,)A, 
(A+ - A_)(4A_A+ + fi,{2 + p,))\h 



ifih - 2A_)(/i, + 2A+) J(2V.(4A_A+ + /i,(2 + p,)) + (1 - A_A+ + ^)A 



^ (A+ - A_)(4A_A+ + /i,(2 + fi,))\f, ■ ^^"^^^ 

The sign will depend on the actual value of the parameters with the requirement that yc 
is real and positive. It follows that in fiuctuations away from 4D Poincare invariance, the 
associated radion has acquired a mass-squared at this order in A. For a large range of 
the parameters, this mass-squared is positive. For example, if fiv,h dominate over A± in 
( p.l6| ), then in this limit the computation of the radion effective potential at second order 
is precisely the one performed by Goldberger and Wise, with the identification of their v^^h 
with our —\v,h/Pv.h- A positive mass-squared results [§|,|lT],|lO . 



The hidden junction condition at y = gives us a fine tuning condition for ph- 
o -o I (2A-(-2 + A_A+)+/i,)A| 

This fine-tuning which we must do is equivalent to fine-tuning the effective four-dimensional 
cosmological constant to zero in order to permit solutions with 4D Poincare invariance. We 
must perform such a fine-tuning order by order in A. 



Let us now check that our basic claims ( p.l7|) are satisfied at this order. For generic 
values of the couplings, 

InS ~ C(l), (2.30) 



so by ( |2.23| ), yc ~ 0{l/e). Given the explicit form for xi it is straightforward to see that the 



resulting A2 satisfies, A' ^ —1. We will show that these successes are maintained at higher 
order in A. 

2. Subdominance of higher orders 



At higher order n > 2, after expanding the square-root in (|2.16|) , we must solve the 
following equation for x 

Xn - 4Xn - eXn " f^vXJ{y - Vc) - f^hXJ{y) = Jn, (2.31) 
where Jn has the form 

•^n = ^ ^r,S,l^ Xh ■ ■ ■ Xi2rXi2r+l ■ ■ ■ X^2{s + r) + l'^*l^ ^-«2{s + r) + l, n- (2.32) 

Here, the expansion coefficients, a are numbers of order one completely determined by 



expanding the square- root in ( p. 16 ), and r + s > 1. Note that J„ is determined from lower 



order solutions of x in perturbation theory. This allows us to treat J„ as a known source in 
the n^^ order equation of motion for x- Thus we can solve ( 2.31|) iteratively for the Xn- 



Our central claim is that to any order in perturbation theory x has the form, 

Xn = E^n>,s,^A"HAz/)"^(Al/c)"^(A/A_)"*e'•^+(^-^=)e^^-^e*^-^^ (2.33) 

where the constant coefficients b are order one or smaller (they can contain positive but not 
negative powers of e), and rii, r, s, t are integers such that 

Hi > 1; r, 712,713,714 > 
71 = rii 7l2 + 713 7I4 

|s|,|t|<7i. (2.34) 
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Let us discuss the significance of tliis claim. First, it shows that the perturbative expan- 
sion is well-behaved. Note that this is a priori not guaranteed. For example, the for- 
mal small parameter A could be overwhelmed by a factor of e^/^ or a high enough power 
of 1/e. However, it is straightforward to see that this is not so given ( p.33| ). Although 
y < He ^ (9(1/A_) ~ 0{l/e), powers of y, yc, and 1/A_ are accompanied by powers of the 
parametrically smaller A. While positive powers of e^+^ > 1 appear, they are compensated 
by powers of 6"'^+^". While arbitrary powers of e^"^ can appear, since we will show (p.l7|) 
self-consistently, e^^'^ ~ C'(l). Secondly, ( p.l7|) also follows since it is clear from ( p.33|) that 



X is dominated by xi for y- Therefore, the second order determination we made for A' 
dominates, and leads to a stable value of yc near that of ( p.27| ). The hidden brane junction 



condition also receives subleading corrections which can be satisfied by fine-tuning ph order 
by order in A. This is the higher-orders incarnation of the four- dimensional cosmo logical 
constant fine-tuning problem of ensuring a 4D Poincare invariant vacuum solution. 

We will now prove our claim ( |2.33| ) for Xn by induction in n. First we note that the 
claim is true for n = 1 based on our explicit solution. For a general perturbative order, n, 
we assume that the claim is true for all lower orders. Then clearly J„, constructed from the 
lower order solution for also has the form 

Jn = Ecn>,s,iA"nAz/)"^(Ai/c)"nA/A_)'^^e''^+(^'-^'^)e^^-^'e*^-^% (2.35) 

where the constant coefficients, c, are also order one or smaller. Given this source term we 
can solve ( p.31|) . 



Xn= r G{y,y')Uy')dy\ (2.36) 

JO 

where G is the Green function satisfying 

1^ - 4^ - e - p,5{y - y,) - /x,5(y)| G(y, y') = 6{y - y'), (2.37) 

subject to orbifold boundary conditions. 

The detailed form of the Green function is straightforwardly worked out. However, in 

order to complete our induction we do not require the full details, but only the general form, 
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o{i)e-^+y' + o{i)e^-(y^-y'^~^+y^j e^+y+ 
o{i)e-^+y' + o{i)e^-^y^-y'^'^+yA e^-y, 



0<y<y' 



oii)e~^-y' + o{i)e-^+y')e^'y, 



+yc)\ p^+y 



(2.38) 



+ 



y' <y < yc- 

The proof that Xn satisfies the claim of ( p. 331) then follows by inspection of all the possible 
terms that can arise from (|2.36|) given (|2.35| ) and (|2.38|) . The basic integrals involved in 
evaluating ( |2.36|) are of the form, 



V2 



dy'y'' 



yi 



or 



/ dyy e y 



yrn+l_yrn+l 

m + 1 



3^5/2 _ p^yi 



(2.39) 



A 



(2.40) 



where m is a non-negative integer, A is some linear combination of A-t with integer coeffi- 
cients, and yi = 0, y; y2 = y, yc- It is straightforward but a little tedious to then check that 
all possible terms arising in ( p.36|) do indeed satisfy the claim of (|2.33|) . 



3. The case of X > e 

We conclude this reworking of the Goldberger-Wise mechanism by discussing the rela- 
tionship between A and e. Thus far, we have considered A, our formal perturbative expansion 
parameter, to be formally smaller than e. The reason for doing this is because of terms in 
( p.33| ) with non-zero n2, or n^, all of which would signify effects of order (A/e)"\ The pres- 
ence of such terms threatens to invalidate perturbation theory if we allowed A > e. However, 
we can in fact prove a stronger result, that all such terms in ( |2.33| ), are in fact accompanied 
by a suppression factor of order e"2+n3+n4^ gQ ^j-^g^^ f^^^j^ take take A > e (though for- 
mally small compared to unity). This relationship is indeed implicit in the original analysis 
of Goldberger and Wise. We have not emphasized this fact until now both for simplicity 
and because the suppression factor is not generally true when we study higher derivative 
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perturbations. We will study the special circumstances under which higher-derivative per- 
turbations do not require A/e to be small in Section V. 

The stronger result, that there are accompanying factors of order g"2+n3+n4^ jg ag^jn 
proven by induction on ra, noting first that it trivially holds for n = 1. For a larger n, we 
assume it holds for lower orders. Therefore, J„ shares the same property since it is made 
from products of (lower-order) x x'- But Xn determined by (|2.36|) can generate (at 
most) one new power of y,yc or 1/A_ relative to J„: a new power of y or yc can arise as 
in (|2.39| ) or (|2.40|) , or one new power of 1/A_ can arise as in (|2.4CI|) in case A = A_. One 
can check (somewhat tediously) that such cases can only arise from terms in J„ involving at 
least one power of x without a derivative acting on it (r > 1 in ( p.32| )). Otherwise, in terms 
where the integrals (|2.39| , |2.40|) produce an extra power of y,yc or 1/A_, the derivatives in 
x' always bring down a power of A_ ~ C^(e) or eliminate a power of y' so that in fact Xn 
does satisfy the claim of the induction. But for terms in J„ with at least one non-derivative 



power of X) we see in ( |2.32| ) that there are explicit powers of e arising in the expansion 
( ^.32| ). These compensate the single power of of 0{l/e) that can be generated in Xn, so the 
induction claim still goes through. 



III. SYMMETRIES OF INTERACTIONS 

Our task now is to specify which higher-derivative interactions are allowed in generalizing 
RSI and the Goldberger-Wise mechanism, in particular how they are constrained by sym- 
metries. Since the spacetime manifold is topologically x S-^ /Z2, the relevant symmetries 
are 5D general coordinate invariance and the Z2 parity symmetry. The coordinate-invariant 
or geometric statement of the Z2 symmetry is as follows. We first start with the manifold 
R^x 5*^ and consider two "3-branes" which divide the 5*^ into two disjoint regions. We choose 
geometries and scalar fields on R^ x such that the two regions are reflection symmetric. 
We now wish to find a convenient description of all interactions which respect this. 

It is useful to start with a formalism which respects full 5D general coordinate invariance 
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on X S^, even though this necessarily involves coordinate systems which do not respect the 



Z2-symmetry of the invariant geometry. (However, see Ref. |T2| for an alternative symmetry 
implementation.) The two fixed-point branes ("hidden" and "visible") will be coordinatized 
as Y;^^{x), Yjll{x), where x are parameterizations of the 3-branes. Using these brane fields, 
a 5D metric, Gmn{X), and 5D Goldberger-Wise scalar x(X), we can form fully 5D general 



coordinate- invariant bulk and brane actions on R^ x S^, as in Ref. |T3 



Once a generally coordinate invariant action has been chosen in this way, we can "gauge- 
fix" by choosing our 5D coordinates, X^'^ , to consist of x^ and an extra-dimensional angle, 
— TT < (f) < TT, such that. 



Y'P — 



yL = ^, (3.1) 
and where the symmetry of the geometry and x is manifest, 

G^,^{x,(\)) = G^^{x, -0) 
Gf,^{x,(j)) = -G^^{x, -0) 
G^^{x, 0) = G^^{x, -(f)) 

Xix,(l)) = x{x,-(f)). (3.2) 

Our procedure for writing actions on R^ x S^/Z2 is therefore to (a) write a general 5D 

coordinate invariant action for bulk and brane fields on R'^ x as detailed in Ref. [113 



(b) gauge-fix the action according to ( p.l|) and (|3.2| ), using the Z2-symmetry of the allowed 
configurations. Once this is done the brane-fields Y,^^{x), Y^l{x), no longer explicitly appear. 
Finally, (c) one must check that all terms in the action are compatible with orbifolding, 
namely they are invariant under 

Gf,^{x,(f)) Gf,^{x, -0) 
Gf,^{x, 0) -G^^{x, -0) 
12 



G^^{x, 0) — > Gcf>^{x, -0) 

x{xA) X{x,-(t)), (3.3) 

for general metrics and scalar field. 

It is straightforward to check that it does not matter whether one imposes the Z2- 
symmetry and gauge-fixing on the action or varies the un-gauge-fixed action without impos- 
ing the Z2-symmetry, as long as one then imposes the Z2 and gauge-fixing on the equations 
of motion. For convenience we will assume that the Z2-symmetry and gauge-fixing have 
been imposed at the level of the action, as implicit in the RSI and Goldberger-Wise papers. 

Recall that the central reason we wish to orbifold is that we wish to take the visible brane 
to have "negative tension" , which would normally yield a ghost-like sign for the associated 
Y^^g kinetic term (Ref. [Q)) causing a vacuum instability towards violent brane fluctuations. 
However, such brane fluctuations violate the Z2-symmetry, and therefore are eliminated by 
orbifolding. It may seem that by formally reintroducing Y^^^ in step (a) of our procedure for 
generating allowed action terms, we are reintroducing the instability. However this is not 
so. The Z2-symmetry renders Y^^^ as pure gauge, with no physical import, as reflected in 
the gauge fixing of step (b). 

IV. HIGHER-DERIVATIVE OPERATORS IN EFFECTIVE BRANE THEORIES 

In any effective field theory we expect there to be higher-derivative operators which are 
remnants of the more fundamental physics which has been integrated out. Their interpreta- 
tion and treatment is generally well-understood. However, they pose special problems when 
they appear in effective theories involving orbifolds and branes in extra dimensions. These 
stem from their appearance as ^-function sources in the extra dimension, rather than some 
structure with a finite size (perhaps due to quantum- mechanical or stringy effects). High 
derivatives applied to such 5-functions produce messy and ill-defined equations of motion. 
We will show how such problems can be treated by a combination of field redefinitions and 
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classical brane renormalization. See Ref. [0 for a recent discussion of brane renormalization 
focusing on co-dimension 2, and Refs. |TH] for earlier work on classical renormalization. 



A. The formal derivative expansion 

To keep only the essential considerations in focus we will not consider gravity here. We 
will limit ourselves to a single scalar field, x? living on an extra-dimensional S'^/Z2 of fixed 
radius (and the usual four dimensions). The same considerations apply in the presence of 
gravity or non-trivial warp-factor. As before we take x to be dimensionless and even under 
the orbifold parity. The general Lagrangian then has the form, 

C = M' l^^idMX? - V{x) - v,{x)S{y-y^) + E ]i^^S.(x,5M)} , (4.1) 

where M is the only explicit scale appearing in order to balance dimensions, and y is now 
defined as the extra-dimensional coordinate corresponding to proper distance in the extra 
dimension. Brane localized terms are represented by multiplying by one of 6{y — yi), where 
yi = 0, 2/2 = ^^c- ^h.d. are arbitrary higher-derivative terms organized in powers of 1/M. 
Note that derivatives always appear in even numbers due to the orbifold parity and 4D 
Lorentz invariance, so that the 1/M-suppressed terms (once the overall is excluded) are 
those with more than two derivatives in the bulk and those with any deriviatves on a brane. 
We have not separated out brane and bulk higher-derivative terms although we will do this 
later. 

We will study the equations of motion subject to the 4D Poincare ansatz, x = xiv)^ 

- dy{K{x)dyx) - -^{dyx? + v\x) + v[{x)S{y -m) = Y. t^^- (4-2) 

This would appear to be a differential equation of arbitrarily high order and therefore requir- 
ing an arbitrary number of initial conditions to solve. However, there is a unique solution 
which is perturbatively close, in a {d/M)"^ expansion, to the zeroth order solution, 

- dy{Kixo)dyXo) - ^^^{dyXo f + V'ixo) + v'Axo)S{y - m) = 0. (4.3) 
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This zeroth order equation is manifestly a second order differential equation subject to 
orbifold boundary conditions, which we know has a distinct solution. Of course, in order to 
determine xo we may have to approximate as we did for the Goldberger-Wise mechanism, 
perturbing in A, but in order to focus on higher derivative perturbations let us take the 
zeroth order solution, Xo, as given. 

An expansion in dy/M may seem sensible in the bulk where solutions are smooth but 
ill-defined in the presence of the brane (5-functions. We will show how these ill-defined brane 
singularities can be renormalized away. To get started and to understand the issues we 
will proceed bravely and formally define the dy/M perturbation expansion.0 We will work 
inductively in n. Let the expansion of x in powers of 1/M be written, 

X = EXm. (4.4) 

m 

Suppose that we can solve ( |4.2| ) to order 1/M™ (we know we can do this for m = 0). We 
will then show how to construct a solution up to order 1/M™+^. 

We substitute ( ^.4| ) into ( [4. 21 ) and focus on the term precisely of order l/M™"*"^, 



{ - dyKixo)dy - dy{dyXo)K'{xo) + lK"ixo)idyXo] 



6S 



+K'{xo){dyXo)dy + V'ixo) + v'lixMy - y^)}xn.+l = ^ 



J2 Xi{y) 



Km 



(4.5) 

m+l 



Sx 

where the right-hand side is the functional derivative of the entire action, including the 
higher derivative terms, evaluated for the field J2e<mXe{y)^ but keeping precisely the terms 
of order 1/M™"^^. Thus, in order to perturbatively improve our solution by one order in 1/M, 
that is solve for Xm+i in terms of Xi<m, we only need to solve the above linear second order 
equation for Xm+i subject to orbifold boundary conditions, with a source term determined 
by the lower order solution XKm- This we can do. 



^This is in the same spirit as setting up the formal Feynman diagram series in quantum field 
theory, which is also ill-defined until after regularization and renormalization. 
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(4.6) 

m+l 

where G is the Green function determined by orbifold boundary conditions satisfying, 

{ - dyK{xo)dy - dy{dyXo)K'{xo) + ^K"{xo){dyXo? 

+K'{xo){dyXo)dy + V'ixo) + v':{xony - y^)}Giy, y') = S{y - y'). (4.7) 

B. The problem 

Such Green functions, G, (and also xo) will clearly be smooth in the bulk, with absolute- 
value type kinks at the branes. That is, their first derivatives will have step-function discon- 
tinuities on the branes and their second derivatives will have 5-functions on the branes. Let 
us suppose this property of the Green function is shared by Xi<m{y) and see how trouble 
can arise in Xm+i- First, consider bulk terms that can appear in ^[Y.e<mXe{y')]\„^+l■ If 
there are more than two derivatives acting on the same field, then we will have derivatives 
of 5-functions on the branes which we must convolute with G, which has absolute- value type 
kinks. The result is therefore ill-defined. Similarly, if we have a power greater than one of 
second derivatives of fields we will have to integrate products of (5-functions which is again 
ill-defined. On the other hand, if there is at most a single field with two derivatives acting on 
it, multiplied by any number of fields with at most first derivatives, then we will only have to 
integrate a single 5-function multiplied by a function with absolute-value type kinks.0 This 
does yield a well-defined Xm+i, which is also smooth except for absolute- value type kinks 
on the branes. It is rather straightforward to see that the troublesome higher-derivative 
bulk terms in the equations of motion correspond precisely to bulk higher-derivative terms 



■^The first derivatives of fields actually have step-function discontinuities, but because of the 
orbifold symmetries they must always appear in even powers, yielding a function with only absolute- 
value type kinks on the branes. 



Xm+iiv) = / dy'G{y,y' 



E xeiy') 

Km 
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in the Lagrangian, ( |4.1| ), with more than one derivative acting on a field (unless it can be 
eliminated by integration by parts). 

Let us now turn to brane terms on the right-hand side of ( |4.6| ), again assuming that 
Xe<m{y) has only absolute- value type kinks and considering the effects on Xm+i{y)- It is 
clear that if in the Lagrangian, ( [4.1D , there are any derivatives in brane terms, then in 
^[J2e<mXeiy')]\,n+i we will have derivatives of 6 functions which makes ([4 .61) ill-defined 
again. On the other hand if there are no derivatives on the brane terms, they will give 
simple (5-functions to be integrated in (|4.6|) , resulting in a Xm+i with only absolute- value 
type kinks. 

Below, we will show that an arbitrary higher- derivative Lagrangian, ( [4.1|) , can be mas- 
saged so that all bulk terms have at most one dy acting on any given field, eg. x^i^xY^, 
and all brane terms have no dy^s at all. As shown above, this will result in a well-defined 
perturbative derivative expansion for solving the equations of motion, with solutions having 
at most absolute-value type kinks. The reader may on a first reading wish to accept that 
this can be done, skip the rest of this section, and continue with the stability analysis for 
higher-derivative bulk terms subject to the above conditions. 



C. Bulk higher-derivative operators 

In this section we will show how higher-derivative operators in the bulk can be massaged 
using field redefinitions. To simplify the considerations note that the equations of motion 
subject to the 4D Poincare ansatz, x = x{y)i can always be obtained by first imposing the 



ansatz on the action, (|4.1| ), and then functionally differentiating with respect to xiv)- Doing 



this we can write the bulk part of the Lagrangian in the form 



C,uik = l-^-K{x){dyx)' - V{x)+J2 

K ft 



^(9xr...(9^xr|, (4.8) 

where \n\ = rii + 2n2... + Nn^ — 2 > 0, and where now d = dy. 

As already discussed, the problematic terms are those with > 1 To begin, we assign 
the bulk potential a formal strength V ~ (9(fc^) < and augment our effective field theory 
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derivative expansion with an expansion in k/M. We then make a field redefinition, ■0(0)) to 
render the kinetic term "canonical" , 

i; = j dxK'/^x)- (4.9) 

Then we have a form 

Uuik = tM' l^(a^)' - V{^) + X: ^(5^)"^-(5^^)"-} , (4.10) 

where of course V and a„ are redefined too. This particular field redefinition will simplify 
stating our procedure, and will be undone at the end. 

We begin by working to zeroth order in k/M (that is, we can neglect V), but to some 
non-zero order, m, in d/M. Working inductively, we assume that by some field redefinitions 
of '0 all > 1 terms have been eliminated from all terms of lower than m-th order. We 
then make the following field redefinition: 

^ - ^ + E (-a)^-^^(a^)-...(a^^)--^ (4.ii) 

\n\=m 

Note that this transformation is only sensible for > 1. To m-th order, the only difference 
this makes to the Lagrangian arises from substituting into the zeroth order kinetic term, 
whereupon (after integrating by parts) it produces precisely the term needed to cancel 
J2\n\=m ^§T^(^'0)"^---(^^'0)"'^ ill the Lagrangian, plus terms of higher than m-th order. Thus 
one uses this procedure to eliminate N > 1 terms to any desired order in d/M, when we 
neglect 0{k'^). 

More generally, the field transformation above will, however, reintroduce terms of various 
orders in d/M, but now with coefficients of order k'^/M^. But we can once again do a field 
redefinition at order k^/M'^ and work to any desired order in d/M to eliminate N > 1 terms. 
This will in turn induce terms of order k'^/M'^. In this way, one can eliminate all > 1 
terms to any fixed order in k'^/M'^ and d/M. Finally, we can transform back to a field, Xnew'- 

^ Ky^iXnew), (4.12) 
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in terms of which, 

ll^m = reM^ \^^K{Xne^){dXne^f - V {Xne^) + E ^]^^(^Xne«,)"} ■ (4-13) 

We have thereby removed all the problematic bulk terms discussed in the previous subsection. 
It can be shown (again, tediously) that our field redefinition in Poincare ansatz, 

x{y)^Xnew{y). (4.14) 

can be lifted to a fully 5D covariant transformation, 

x{^,y)^Xnew{x,y), (4.15) 

which, however, reduces to x(l/) ~^ Xnewiu) upon imposing the 4D Poincare ansatz. Simi- 
larly, when gravity is included there are fully 5D generally covariant field redefinitions which 
upon imposing the Poincare ansatz ensure that all > 1 terms (terms with more than one 
dy acting on a field) are absent. 



D. Brane higher-derivative operators 

It is obvious that the field transformations which we discussed above to massage the bulk 
action will induce derivative terms on the orbifold fixed points. Furthermore, such derivative 
terms may already be present anyway. Therefore after ensuring that all bulk terms satisfy 
A^" < 1, we must still massage away brane terms in the Lagrangian of the form, 

Uane = / dyS{y - y,)M'k i^ ^. ^(dxr-id^xT^^ ■ (4-16) 

We will show that the important physics contained in x away from the branes can be 
obtained by using effective renormalized brane Lagrangians without any extra-dimensional 
derivatives, 

jOeff. brane = J dy6{y - yi)M^kvi{x). (4.17) 
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The first step before renormalizing tlie original Lagrangian is to regulate it, say by the 
replacement, 

Ky - Vi) D{y - = Ne-^^^y-y^+'^e^/^y-y^-'\ -t<y<t 

0, else, (4.18) 

where 

iV-i = rfye-i/(-^+*)ei/(^-*). (4.19) 

This replaces the infinitesimally thin brane by a thick brane of thickness t: 1/M < t << r^. 
Note that the brane profile D{y — yi) is a smooth function (that is, a C°° function) of 
compact support. 

With this regulator in place, clearly our perturbative derivative expansion is always well- 
defined despite derivatives on branes, and ( [4.6[ ) will always give smooth solutions. Of course, 
the results depend on our choice of regulator. However we will show that away from the 
(compact) core of the branes, this regulator-dependence can be completely subsumed into 
an effective brane Lagrangian, ( ^.17] ). 



For simplicity of exposition we will consider a single orbifold fixed point on R/Z2 rather 
than the two fixed points of jTji- The generalization to 8^/7^2 is entirely straightforward. 
First, recall that we have shown that at every order in perturbation theory we only solve 
second order differential equations. Thus the solution xiv) perturbation theory is com- 
pletely specified by x(0) and dyxiS^)- Orbifold parity sets 9yx(0) = 0, but x(0) (on R/Z2 
not S'^ /7i2) is an unfixed number, x(0) = c. In particular by solving the equations of motion, 
x(t) = x(— t) and dyxit) = —dyx{—t) are both functions of c, 

X{t) = fic) 

dyx{t)=g{c). (4.20) 

Let us define xiv): to be the result of integrating the second order differential equations from 
t to a general y > 0, using xit)^dyx(t) as boundary conditions, but completely neglecting 

20 



brane interactions. Now for y > t, x{y) is in fact the correct solution x(y), since the 
neglected brane interactions vanish in this region. But for y < t, clearly xiu) cannot be 
trusted. Nevertheless, as long as the physics of interest is dominated by bulk field behavior 
outside the core of the brane, we can use x(y). This will be the case for the solution to the 
hierarchy problem, which is determined by the warp factor accumulated over the large bulk 
(re> 1/M). 
Let us define, 

dyx{0+)=G{x{t),dyx{t)). (4.21) 

We can eliminate dependence on x(i), dyx{t) in favor of c, 

x(0+) = F(/(c),^(c))=p(c) 
a,x(0+) = G(/(c),^(c))^g(c). (4.22) 

Therefore, order by order in perturbation theory we can eliminate dependence on c by 
inverting p, 

5.X(0+) = g(p-'(x(0+))). (4.23) 

We will now show that x is the classical solution corresponding to a "renormahzed" ef- 
fective Lagrangian with the same bulk terms as before, but with a 5-function brane- localized 
potential term (that is without any extra-dimensional derivatives). Recall that by field re- 
definitions we have already ensured that the bulk Lagrangian only depends on x and its first 
derivative, >Cb„;fe(x, x')- The full effective Lagrangian, including the effective brane potential, 
is then given by, 

Ceff = C^lkiX. X') + <^(z/)^e//(x), (4.24) 

where, 

Veffix) = jdx MP-\X))^^ {X,q{p-\X))) . (4.25) 
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The reader can straightforwardly check that the equation of motion that follows from this 
effective Lagrangian is equivalent to the equation of motion due only to the bulk terms away 
from y = 0, supplemented by the boundary condition, ( [4.23D . 

Thus, we are always able to find an effective Lagrangian of the form ( |4.24| ), which has 
solutions which agree with those of a general Lagrangian supplemented by a regulator, 
outside the regulated core of the branes. The important physics such as the hierarchy are 
insensitive to the general disagreement inside the thickness of the regulated brane. Note 
that the effective Lagrangian does not suffer from UV ambiguities, and therefore needs no 
regulator. In that sense it is "renormalized" . 



V. STABILITY AND HIGHER DERIVATIVE OPERATORS 
A. Set-up of the model 

We will now show that the Goldberger-Wise mechanism can be realized in a very general 
setting, including higher-derivative operators. We will take the higher-derivative terms to be 
suppressed by appropriate powers of 1/M and constrained by symmetries according to the 
discussion of Section III. Furthermore, we will assume that they have already been massaged 
by field redefinitions and brane-renormalization as discussed in the previous section, so that 
fields appearing in bulk interactions have at most one extra-dimensional derivative acting 
on them, while there are no extra-dimensional derivatives in brane terms. 

The action in Einstein frame is 

Sbulk = M' J d'x jy^pVG (^eV{x) - \R+lK{x){Dxy + M'C,.,XGmn,X,^)^ (5-1) 

S,is = -M^k J (fx^sv.ix) (5.2) 

Shid = -M^k J d'^x^/gj~ivhix)- (5-3) 

Here, Ch.d. indicates terms containing more than two derivatives with dimensionless coeffi- 
cients of order unity and suppressed by powers of 1/M. We can expand the dimensionless 
functions V, K, v^, and Vh in powers of x- 
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V{x) = 3 - \tx^ + 0{x') 
K{x) = 1 + Oix) 

Mx) = Pv + Kx + \pvX^ + 0{x^) 

Vh{x)= Ph + Kx + \phX^ + 0{x^). (5.4) 

We will perform perturbation theory in powers of k^/M^ <^ 1 and the brane-tadpoles, 
Xy^h, as discussed in Section II. Higher-derivative terms are automatically suppressed by 
powers of k'^/M'^ because the bulk potential is dominated by a 5D cosmological constant of 
order M^k'^. Formally, this is seen by the fact that in our dimensionless y-coordinate every 
dy is accompanied by k. For simplicity, we take k'^/M'^ and X^^h to all be of the order of 
our formal small parameter, A. We continue to take the relations of ( |2.6| ) to hold among 
Xv,h) Pv,hi Pv,h and e. 

We now write the equations of motion with the Poincare ansatz using again the dimen- 
sionless variable y = krc(p as in Section II, 

X" + 4AY - ex = (A, + + Oix'mv - yc) + (A^ + PhX + 0{x^my) 

+ f,{x,x\A') + x"f2{x,x\A') (5.5) 

A" + lixT = -liPv + A.X + Ipvx' + Oix'My - yc) - ^(p. + A,x + ^Phx' + Oix'My) 
+ Ux,x\A') + A"U{x,x',A') (5.6) 

{AT = 1 - ]:^X' + kxT + Hx, X, A'), (5.7) 
6 6 

where the refer to the variation of the bulk higher- derivative action subject to the special 



form discussed in Section IV. Note that the right-hand side of (5]7) does not contain terms 
with delta functions or second derivatives. 

Perturbation theory follows along lines similar to Section II. We first formally expand A 
and X as in ( p^.l4 ) and substitute into the equations of motion. We then solve ( p.7[ ) for the 



A'j^ in terms of the Xm<n and x'm<n to any desired order. Substituting for the A'^ into 
then yields an equation involving only the Xn and x'n which we solve perturbatively. Then, 
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by 5D general covariance in the bulk, ( p.6|) is automatically solved, except for the two brane 
junction conditions which are solved by fine-tuning and setting Hc to its stable vacuum 
value. 

To see that ( |5.6| ) is automatically satisfied up to junction conditions, note that under 
infinitesimal general coordinate transformations the bulk action is invariant. 



bulk 



Gmn + Dmtin + Dntim, X + d'^XVM = Sbuik [Gmn, x] , (5-8) 



for arbitrary infinitesimal %/. This implies 



In the Poincare ansatz this gives 

^ I SSbulk\ ./SSbulk 1 i^Sbulk /r in\ 

This shows that (|5.6| ), the G^u equation of motion, is satisfied up to junction conditions, 
given (|5T5|), the x equation of motion, and (p.7|), the G55 equation of motion. 



The equation for x obtained by eliminating A' from ( |5.5| ) using ( |5.7] ) has the form ( p. 311 ), 



but now with sources following from ( p.5[ - p.7| ) of the form 



where 



Jn = Jn + Jn + (5- H) 

Jn'" = E %s,^^"'X^, • • • XvXl,, • • • x'.A+-+^s+m, n (5.12) 
jtrane, ^ ^ 6 .A^Xn • • • X^A.+ ^.+m,nS{y - %), (5.13) 



and 



jmixed \ ^ 



r,s,?'^ ^ii ■ • • XiaXi^^i ■ ■ ■ XigXp^ii+---is+m+p, n- (5-14) 

This last term acts as a brane term as well as a bulk term because x" contains (5-functions 
at the branes as well as a smooth bulk behavior. Again, the constant coefficients, a, h and 
c are order one or smaller. The n*^ order solution is obtained using ( f^.36D with the same 
Green function defined by (|2.371 ). 
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B. Stabilization at second order 



We will show here that hierarchy stabilization satisfying ( ^.17] ) is naturally achieved by 
second order in A, similarly to Section II. Note that at zeroth order our solution is again 
the unstabilized RSI solution. At first order xi is given again by ( |2.19| ), although A[ may 
now be a non-zero constant. As discussed above, we can solve for the A'^^ using (|5.7|) in 
terms of the Xm<n and Xm<n- Therefore, A2 must be a quadratic polynomial in xi and Xi- 
Summarizing, we have 

A'q + A[ + ^2 = 00 + ttiXi + «2xi + ^ex? - ^xf , (5.15) 

where ao = —1 + 0{X), a\^2 ~ ^i}^- The constants ai are independent of yc since the 
are determined in terms of the Xm. and Xm locally. 

It remains to satisfy the junction conditions for ( |5.6|) . After an 0(A^) tuning of (not 
appearing at higher orders), similar to Section II, the above considerations and ( p.l9|) yield 
a visible junction condition of the form 

P2(e^-^=) = 0, (5.16) 

where P2 is a quadratic polynomial with order one coefficients. For generic values of these 
coefficients one obtains solutions, 

e^-y- = 0{1). (5.17) 

Thus, along with ( p.l5| ) we have demonstrated the hierarchy, (p.l7|) , at this order. The 
hidden brane junction condition is solved by fine-tuning ph, corresponding to the fine-tuning 
of the 4D cosmological constant to zero. 

C. Subdominance of higher orders 

We claim that at any order Xn has the general form ( p.33| ) , and A'^ has the general form. 



K = E9n>,s,tA"HAy)"^(Aye)"^(A/A_)"^e'^^+(^-^^)e^^-^e*^-^^ (5.18) 



n / y 
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so that perturbation theory is well-behaved. Therefore, our second order mechanism above 
for stabilization of the weak/Planck hierarchy of order e"*^*^^-*/^ continues to hold. 

We show this claim is true by induction. We have already discussed the cases for n = 0, 1. 
Focusing now on Xn and assuming that ( |2.33| ) is valid at all orders lower than n, the n^^ 
order source has the form 

+ E^n>,«,t^"HAl/)"^(Ayc)"^(A/A_)"^e-^^+^^e*^-2'=5(l/) 

+ Ecn>,.,tA"HAl/)"^(Ai/e)"^(A/A_)"*e-^^+^=e*^-^=5(2/ - y,). (5.19) 

It is straightforward to check that with this source, the solution for Xn given by ( p.36[ ) indeed 
satisfies ( p.33|) . Then, given the form of (|577|) and the form for Xn of (|2.33| ), ( ^.18|) readily 



follows. 



D. The case of A > e 



A significant difference between the source considered in this section, ( |5.11|) and that 



considered in Section II, ( |2.32|) , is that in the latter every power of x without a derivative 
was accompanied by an explicit power of -y/e. Recall that this ensured that increasing powers 
of 1/A_ in perturbation theory were accompanied by e factors. This allowed us to take A > e 
and still have a meaningful perturbative expansion. Without this feature, in this section we 
are limited to formally A < e, as we have assumed thus far in this section, in order to have 
a good expansion. 

However, there is an interesting scenario under which we would recover a good expansion 
for A > e even for higher- derivative operators. We first assume a (non-linearly realized) 
global symmetry^, 

X — > X -|- constant. (5.20) 



^We thank W. Goldberger for pointing out the usefulness of this symmetry to us. 
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If this symmetry were exact, only derivatives of x would be permitted in the action. But 
one does not expect that such global symmetries are respected by the underlying quantum 
gravity theory We will assume that in fact such violation is controlled by the small param- 
eter, e, so that non-derivative powers of x in the action are accompanied by \fe^ thereby 
generalizing the explicit appearance of e in Section II. Then, increasing powers of 1/A_ in 
perturbation theory are again accompanied by e factors, so that the perturbative expansion 
is under control for A > e. 

VI. CONCLUSIONS 

We demonstrated the stability of the RS mechanism for generating the weak/Planck 
hierarchy in the presence of higher-derivative interactions. We re-worked the Goldberger- 
Wise radius stabilization mechanism in a systematic perturbative expansion in parameters 
of the brane potential. Incorporating higher- derivative interactions as further perturbations, 
we showed that they did not affect the basic mechanism. 

In our perturbative analysis, radius stabilization is achieved in the following steps. At 
zeroth order we found the unstabilized RSI vacuum with a trivial profile for the Goldberger- 
Wise scalar. At first order, the scalar responds non-trivially to the brane-potentials. The 
back-reaction of this scalar profile on the 5D metric (warp factor) takes place at second 
order. In particular, the radius acquires a distinct, stable value. This value can naturally 
be several times the fundamental scale of the theory, while the bulk geometry is a small 
perturbation of AdSs, so that the RS mechanism for generating the hierarchy operates. We 
also gave a careful treatment of all higher orders in perturbation theory, showing that our 
expansion is under good control. 

The derivative expansion of effective field theory seems at first fundamentally at odds 
with the presence of "thin" or 5-function branes or orbifold fixed points, resulting in ill- 
defined singularities in the classical equations of motion. While such singularities can be 
regulated, the process is messy and apparently regulator-dependent. However, we showed 
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how classical renormalization of the brane action can be implemented so as to remove the 
need for explicit regulators. This greatly simplifies the discussion of higher-derivative effects. 

It is important to study the stability of the RSI effective field theory under quantum 
effects. Since the effective theory is not renormalizable we expect quantum divergences of 
the form of every possible local operator (subject only to symmetries). However, once these 
divergences are covariantly regulated at the fundamental scale, they must be of the form 
and strength of the higher- derivative interactions already considered in this paper. As usual 
in effective field theory, renormalization proceeds order by order in the derivative expansion, 
using counterterms also of the form of the operators of this paper. Thus, because quantum 
divergences are local they cannot destabilize the RS hierarchy, since our purely classical 
analysis already treats all such local effects. This leaves only the non-local quantum effects 
which are UV-finite and therefore well-defined and calculable. In future work, we hope to 
build on Refs. in studying the general structure of these non-local quantum amplitudes 
and whether they affect the basic RSI mechanism. 
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